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We evaluate the electronic transmission and conductance in bilayer graphene through a finite 
number of potential barriers. Further, we evaluate the dispersion relation in a bilayer graphene 
superlattice with a periodic potential applied to both layers. As model we use the massless Dirac- 
Weyl equation in the continuum model. For zero bias the dispersion relation shows a finite gap 
for carriers with zero momentum in the direction parallel to the barriers. This is in contrast to 
single-layer graphene where no such gap was found. A gap also appears for a finite bias. Numerical 
results for the energy spectrum, conductance, and the density of states are presented and contrasted 
with those pertaining to single-layer graphene. 
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I. INTRODUCTION 

Low-dimensional systems have long been the subject 
of intensive research, both on their fundamental prop- 
erties and on possible applications. In this respect the 
recent production of atom-thick crystal carbon layers 
(graphene) has raised the possibility of the development 
of new graphene-based devices that exploit its unusual 
electronic and mechanical properties (for a recent review 
see Ref . [T| . The electronic spectrum of defect-free single- 
layer graphene is gapless and, together with the chiral 
aspect of the carriers in this system, leads to a perfect 
transmission through an arbitrarily h igh and wide po- 
tential barrier, i.e., the Klein paradoxPS. That can be 
avoided if a gap is introduced in the electronic spectrum 
and may be necessary for certain applications, e. g., for 
improving the on/off ratio in carbon-based transistors. 

There are a few methods to introduce a gap in the 
spectrum of graphene. One of them is to use nanoribbons 
in which a bandgajP arises due to the lateral confinement. 
Also, depositing graphene on a substrate such as boron 
nitride was found recently to result in a bandgapP of 53 
meV. In bilayer graphen^ a gap can be introduced by 
applying a bias between the two layers or by doping one of 
them such that a potential difference results between the 
layer^^KlS]^ Changing the bias in the latter case can open 
and close the gap dynamically which is interesting for 
transistor applications. Nanostructured gates can thus 
allow the creation of quantum dots on bilayer graphen^^^'. 

In this paper we investigate the electronic properties of 
a biased bilayer in which the potential difference between 
the two layers is changed periodically. Such a superlat- 
tice (SL), which can be created by applying gates to the 
bilayer, is of interest as it shows how a one-dimensional 
band structure may appear in such a system. An addi- 
tional motivation is that curvature effects of corrugated 
single-layer graphene lead to an effective periodic poten- 
tial resembling that of a SLp^. Although in a bilayer this 
effect would be weaker, since the bilayer is less bendable 



than a single layer, it might still be important. 

The paper is organized as follows. Section II briefly 
shows the basic formalism. In Sec. Ill results for the 
transmission and conductance through a finite number 
of barriers are presented. Section IV shows results for 
the dispersion relation and the density of states in SLs 
in bilayer graphene. Finally, a summary and concluding 
remarks are given in Sec. V. 



II. HAMILTONIAN, ENERGY SPECTRUM, 
AND EIGENSTATES 



Bilayer graphene consists of two A-B-stacked mono- 
layers of graphene. Each monolayer has two independent 
atoms A and B in its unit cell. The relevant Hamiltonian, 
obtained by a nearest-neighbour, tight-binding approxi- 
mation near the K-point and the eigenstates 4* read 
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Here tt = {p.j., + ipy), 



-ihdx y is the momentum 



operator, vp ~ 10 m/s is the Fermi velocity, Vi and V2 
are the potentials on layers 1 and 2, respectively, and tj_ 
describes the coupling between these layers. As shown 
in Appendix A, for spatially independent t±, Vi, and V2 , 
the spectrum consists of four bands given by 
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Here e^^, — P 



+ A = (Vi - V2), 5 = A/2hvp, 
£ — E/hvp and t' — t±/hvp. The eigenstates ip oiH are 
given by Eq. ( 32 1 in the Appendix. 
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FIG. 1: (Color online) Transmission through a single square 
barrier of height V = 100 meV and width 10 nm. The po- 
tential difference (5 = (Vi — V2)/2ft«F is zero in panel (a) and 
A = 100 meV in panel (b), inside the barrier/well region. 
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FIG. 2: (Color online) Transmission through a double bar- 
rier. The square barriers are 100 meV high and 10 nm wide, 
the distance between them is L = 10 nm. In panel (a) the 
potential difference between the layers is zero, in panel (b) it 
is 100 meV inside the barrier /well regions. 



potential barrier. The eigenstates ^ given in Appendix 
A can be used in each region of constant potential. In 
matrix notation the wave function in region j with a con- 
stant potential can be written as a matrix product (cf. 
Eq. (l34|) of the Appendix), 



(5) 



where Aj — [Aj, Bj, Cj, Dj]^ and the superscript T de- 
notes the transpose of the row vector. Then we apply the 
continuity of the wave function at the different potential 
steps. For the {j + 1)*'* potential step at Xj+i we obtain 



(6) 



This links the coefficients of the wave function behind 
the barriers to those in front of them. Then we can write 
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where Afj = Mj^i{xj+i)Qj'_^iQjMj{xj+i) and Af = 
Yij-^j- From now on we assume \E\ < t± outside the 
barrier such that a+ S K and a_ G C, see the Appendix. 
Assuming that there is an incident wave, with wave vec- 
tor from the left (normalized to unity) , part of it will 
be reflected (coefficient r) and part of it will be transmit- 
ted (coefficient t). Also there are growing and decaying 
evanescent states near the barrier (coefficients eg and e^, 
respectively). The relation between all these waves is 
written in the form 
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A reduced version of the four-band Hamiltonian shown 
in Eq. (1) that is often usecP is given by 
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Assuming solutions of the form Aeiqi{ikxx) exjp{ikyy) we 
can replace Px by hk^- Then setting the determinant of 
the equation Ji^p = Eip equal to zero gives rise to the 
two-band spectrum 



E-V = ±{vF^h^/t^){kl + kl), 



(4) 



where V is the potential applied to each layer. In the 
next section we compare some of the results obtained 
from Eqs. (1) and (2) with those obtained from Eqs. (3) 
and (4). 



III. FINITE NUMBER OF BARRIERS 

A. Transmission 

When applying a bias to a metallic strip a potential 
barrier is created; we will approximate it by a square 



it can be rewritten as a linear system of equations, 
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where Nij are the coefficient of Af. We solved this set of 
equations numerically. The transmission is now given by 
T=\t\\ 

A contour plot of the transmission through a single 
barrier is shown in Fig. [l] Panel (a) is for a barrier with 
height 100 meV and width 10 nm and the potential dif- 
ference A between the layers is zero. In contrast with 
the case of a 2DEG, there are transmission resonances 
in the region which corresponds to energies lower than 
the barrier height. These are due to the hole states in- 
side the barrier through which the electrons can tunnel. 
In panel (b) the barrier is 10 nm wide and the potential 
difference between the layers is A = 100 meV. As can 
be seen, the ky dependence of the transmission in panel 
(b) is weaker than that on k^ and resembles more the 
Schrodinger case. The transmission of electrons, at nor- 
mal incidence {ky = 0), starts from a, kx — 0.23 nm~ 
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FIG. 3: (Color online) Transmission through a 100 nm wide 
barrier as a function of the angle of incidence for constant 
energy E = 17 meV. Panel (a) results from the 4x4 Hamil- 
tonian of Eq. (1) and panel (b) from the 2x2 one of Eq. (3). 
The solid red and dashed green curves are for a single barrier 
with height 50 meV and 100 meV, respectively. 
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FIG. 4: (Color online) Conductance as a function of energy. 
The very thick blue, thick red, and thin black solid curves are 
for two, five, and ten barriers, respectively, of width 10 nm, 
height 100 meV and with an interbarrier distance of 5 nm, 
and result from the 4x4 Hamiltonian, Eq. Jll, while the 
dashed curves result from the 2x2 one, Eq. (IsTT 



B. Conductance 



meV which corresponds to an energy of 50 meV which 
is the edge of the gap, inside which there are evanescent 
states which supress the transmission. For the double 
barrier system we see that there is also a resonance at 
kx — 0.16 nm~. In contrast with the single-layer case, 
in Fig. 1 (a) and 2 (a) there is no perfect transmission 
for normal incidence {ky = 0), even though the system is 
gapless. This is a consequence of the chiral nature of the 
carriers in bilayer graphene (see, e.g. Ref. 

A plot of the transmission through a double barrier 
is shown in Fig. [2] The barriers are 100 meV high, 10 
nm wide, and the distance between them is L = 10 nm. 
Panel (a) is ioi S — and panel (b) for A = 100 meV. 
In agreement with Ref. 13 we find that it is the distance 
L between the barriers and not their width that is im- 
portant in determining the tunneling states and thus the 
transmission. For more results, e. g. conduction through 
unbiased, multiple-barrier systems see Refs. [T51and UM 

The transmission shown in Fig. [T] and [2] depends on 
the angle of incidence (j) given by tan0 = ky/kx- A more 
direct way to see this is shown in Fig. |3] where the trans- 
mission is plotted as a function of the angle of incidence 
for constant energy E — YI meV. For panel (a) we used 
the 4x4 Hamiltonian of Eq. (1), where a bias A = 0, 
while for panel (b) we used the 2x2 one of Eq. (3). The 
differences are due to the different spectra of Eqs. (2) 
and (4). We consider the results of panel (a) to be more 
accurate than those of panel (b), since the Hamiltonian 
of Eq. (1) includes the effect of the non-parabolicity of 
the electron dispersion, which can have a strong effect on 
the resonant transmission of the carriers. 



It is interesting to see to what extent the transmission 
affects the conductance G, which is given by 

G = Go T(£;,(/))cos0#. (10) 

J-Tr/2 

Here Gq = 2e^ -I- t^^EpLy/ {-nkhvp), 4> is the angle 
of incidence measured from the x axis, T{E, </>) the trans- 
mission through the structure at energy E, and Ly the 
length of the structure along the y direction. 

In Fig. |4] we plot the conductance G through two, 
five, and ten barriers in blue, red, and black color, re- 
spectively. The height of the barriers is 100 meV, their 
width D = 10 nm, and the interbarrier distance L = 5 
nm. The solid curves are obtained using Eq. (1) and the 
dashed ones using the reduced Hamiltonian of Eq. (3) 
with the same coupling strength t± = 390 meV. As can 
be seen, both models give qualitatively the same results. 
The disagreement is mostly apparent in the low-energy 
region and is mainly due to the large deviation E — V 
of the energy from the barrier potential V , due to which 
the 2x2 Hamiltonian approximation inside the barrier 
fails. 

IV. SUPERLATTICE 

A. Dispersion relation 

The model we used for a superlattice (SL) in graphene 
is shown schematically in Fig. [5] The electronic spec- 
trum resulting from this periodic structure can be ob- 
tained by writing the solution for the spinors as Bloch 
waves and applying the continuity condition for the wave 
function at the potential steps. In both barriers and 
wells the solutions are the ones for a constant poten- 
tial and the boundary conditions determine the matrix 
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FIG. 5: (Color online) Schematics of two experimental setups 
for realizing the three SL potentials we investigated. In panel 

(a) the layer potentials Vi and V2 are kept the same, the 
experimental setup shown can be used. The setup in panel 

(b) can establish a bias A — Vi — ¥2- In both experimental 
setups the layer potentials are controlled by the applied top 
Vtg and back Vtg gates. 



relation between the wavefunction coefRcients in the two 
regions. For a periodic potential, Bloch's theorem applies 
with period I ~ a + b, implying ipk{x + I) = 4'k{x)e^^^ . 
Then referring to Fig. [5] we obtain: ^"(0— ) = ^^(0+) 
and tp{—a) = '0(6)e~*'^'. Writing the wave function in 
the regions of constant potential as a matrix product 
\]/ = QAiA, labelling the coefficient matrices inside the 
barrier regions as Ai and the ones inside the well regions 
as A2, and applying the above boundary conditions we 
obtain the matrix equations 



GiAi = g2A2 

giMi{-a)Ai = ^2X2(6)6-^^^2. 
Eliminating Ai in Eqs. (10) and (11) leads to 



Equating the determinant of Eq. ( 13 I to zero 
det[Xi(-a)gfig2 -ar'e2A^2(6)e-'"] = 



0. 



(11) 



(12) 



(13) 



(14) 



The solution of Eq. (14) gives the energy spectrum or 
dispersion relation. From this determinant we search for 
the zeros of Eq. (14) using the Newton method and obtain 
the dispersion relation. 

In Fig. [6] we plot the dispersion relation versus and 
ky for three different SLs. In the first one we take the 



potential on the back and front gates to be the same V = 
50 meV; between the strips the potential is —50 meV. In 
the other two we only vary the bias difference A between 
the two layers: A = 50 meV in the barriers and meV in 
the wells for the second SL and correspondingly A = 50 
meV and 25 meV for the third one. The average potential 
of both layers is kept constant. The parameters used are 
a = 5 = 10 nm and the tunnel coupling t± = 390 meV. 
Only the first two minibands are shown in the left panels 
of Fig. 6. The middle column shows cross sections of 
the dispersion relation for constant kr^ in panels (a) and 
constant ky in panels (b). 

Applying the first type of SL potential shows the for- 
mation of subbands. For this first type one can easily find 
an analytical expression for the one-dimensional case, i.e. 
ky = 0, the formula for the dispersion calculated from the 
four-band Hamiltonian is (union of formula with + and 



cos(fcL) = cos(a± ifl) cos(q;± 26)— G± sin(Q;± la) sin(Q;± 2b) 

(15) 



Where £j = E — Vj in region j 
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The analog of this formula for the two-band approx- 
imation of the Hamiltonian is the same formula with 
a±j = -^iej instead of a± j- = -^/e^ ± £jt±_ 

In contrast with the gapless spectrum of SLs on single- 
layer graphene, here a bandgap is found for ky — 0. This 
is in agreement with the fact that for the transmission 
through a barrier there is no perfect transmission for per- 
pendicular incidence while in single-layer graphene there 
is. The second SL potential has the same barrier and 
well parameters as Fig. 2(b), the resonance at kx — 0.16 
nm^^ meV which we saw in this double barrier system 
seems to correspond to the energy value of the first band, 
E ~2h meV. Also, for the lowest band the mexican-hat 
energy profile of biased bilayer graphene is retained in the 
ky direction. In the third SL potential the gaps between 
the subbands are smaller than those of the second SL and 
the dispersion relation resembles more the (folded) one of 
bilayer graphene without any SL potential but with an 
applied constant potential difference. The DOS of the 
latter two SLs shows large van Hove peaks at energies 
corresponding to the lowest band, there also the velocity 
is zero and localized states form. 

B. Density of states 

To understand part of the behavior of carriers in a SL 
we evaluate the density of states (DOS) D{E). In the 
reduced-zone scheme it is given by 



TV/I 



dkySiE~E„ikx,ky)), (17) 



5 



A V,(x) = V,(x) 

jTuin. 




ijj-lIj-LJ"- 



> 
E, 

LU 



-75 
0.2 



0.1 
k (nm"' 



0,5 

k (ti/I) 



^^^^ 













0.2 0.1 0.0 0.5 1.0 

k(nm-') k(„/l) 



0.2 0.1 0.0 0.5 1.0 

k (ran"') k (nm"') 




-200 -100 100 
E (meV) 



-300 -200 -100 



100 200 300 



E (meV) 




-300 -200 -100 100 200 300 

E (meV) 



FIG. 6: (Color online) Dispersion relation and DOS for three types of SLs. (1) The barriers are 50 meV high and the wells —50 
meV deep. (2) The barriers are biased by A = 50 meV, the wells are unbiased, i.e., A = meV. (3) The barriers are biased 
by A = 50 meV and the wells by A = 25 meV. Left column: energy vs fe^, and ky for a = 6 = 10 nm and tx_ = 390 meV. 
Lines of constant energy, belonging to the lower miniband, are projected onto the {kx,ky) plane. Middle column: slices of 
the corresponding dispersion relation, (a) for constant k^ = (solid magenta curves) and k^ — n/l (dashed green curves), and 
(b) for constant ky = (solid red curves) and ky = 0.2/nm (dashed blue curves). Only half the Brillouin zone is shown. Right 
column: DOS for the corresponding SL. For the unbiased SL (1) we also show the DOS (red area) for a SL with the same 
parameters on a single-layer graphene. The dashed (dash-dotted) curves show the bilayer (single- layer) DOS in the absence of 
the SL potential. 



where A is the surface area. The integral is evaluated where the and ky indices take the values 
numerically by converting it to a sum in the manner 



/ dkx I dky ^ [ ) ( """^ ) \^ , (18) The cutoff fcmoK for fcy is chosen sufSciently large, we took 

'0 •'0 Ny J. _g J, _Q fc™„„ = 2 nm"-"^. Tn addition, we renlace the S function 
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in Eq. (15) by a gaussian, 

S{E-E^{kx,ky)) « (1/V2^) e-^E-E„iK,k,r/2a^^ (20) 

and choose a small but sufficiently large to compensate 
for the discretization of k^ and ky, i. e. we took cr = 0.03 
mcV. The evaluated DOS is shown in the right panel of 
Fig. 6. In these figures the magenta, green and orange 
areas are for bilayer SLs and the red one for a single- 
layer SL. The dashed and dash-dotted curves show the 
DOS for single-layer (Ds) and bilayer (Dt,) graphene in 
the absence of the SL potential given by 



D,{E) = \E\/hvF, 

DbiE) = i\E\+t^/2)/hvF, 



where we used the usual tight-binding HamiltoniarP^ for 
single-layer graphene and the one given by Eq. ([T]) for 
bilayer graphene. The peaks in the DOS have the typical 
1/^/E — Eq behavior of ID subbands. 



V. SUMMARY AND CONCLUDING REMARKS 

We evaluated the electronic transmission and conduc- 
tance through a finite number of bilayer graphene barri- 
ers. Further, we obtained the dispersion relation and the 
DOS for a periodically biased bilayer, i.e., a bilayer in 
the presence of a SL potential. With the rapid progress 
in the field we expect that such a periodic biasing will 
soon be realized experimentally. Since the elastic mean 
free path of carriers in high-mobility graphene layers can 
be of the order of hundreds of nanometers, a ballistic be- 
havior can be expected to be observable on the length 
scale of the periodic structures discussed here. 

For some transmission and conductance results we used 
both the four-band Hamiltonian given by Eq. ([T]) as well 
as the reduced two-band Hamiltonian given by Eq. (3), 
cf. Fig. 4 and 5. We consider the former results as 
more accurate than the latter ones, since the graphene 
bilayer spectrum obtained from the four-band Hamilto- 
nian is known to give a better agreement with both exper- 
imental data and theoretical tight-binding calculations^. 

For zero bias the dispersion relation shows a finite gap 
for carriers with zero momentum in the direction par allel 
to the barriers in contrast to the well-known resultJ^^ 
for single-layer graphene, cf. Fig. 6. A gap also appears 
for a finite bias, cf. Fig. 6. We also contrasted the 
DOS for bilayer graphene with the corresponding one for 
single-layer graphene, cf. Fig. 6. We expect that all these 
results will be tested experimentally in the near future. 
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VI. APPENDIX 

We assume solutions of the form 'ifc{x,y) = 
(j)c[x)e'''yy,C = A,B. Then Eq. (1) and Schrodinger's 
equation Hil) = E^j lead to the following equations 

-i{dx - ky)(j)B - (e' - S)(I)A - t'(t)B', (22a) 

-i{dx + ky)^A = {e' -S)(I)B, (22b) 

^t{dx + ky)(j)A' - (e' + S)(bB' - t'^A, (22c) 

^i{dx-ky)^B' ^{e' + S)(bA', (22d) 

where e' ^ 6 — e — {uq ± (5), ui — uq + S, and U2 — uq ~ S. 
We solve Eq. (|22b|) for 0b and Eq. (|22d|) for and 



(21) substitute the results in Eqs. (22a) and (22c). This gives 



(dl - kDct^A = -{e' - 5fct>A + t'{e' - S)q^, 



(23a) 



{dl - kl)(t>B' = -{e' + 5f<l)B' + t'{e' + 6)(I)A- (23b) 

For the system of Eqs. ( |22| ) and for constant potentials 
the spectrum is determined by the equation 

[-k^ + {e'-Sf][-k'' + {e' + Sf]-t'^{e'^-S'') = 0. (24) 

Solving it leads to four bands (e^j, = k^ + + 
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(26) 



tit, ± iV4fc2^Vi'2 + fc2 + i'2/4 
and four possible wave vectors ±a± — (fc^ — ^y)^^^ 

1/2 



(27) 



To obtain the general solution for the spinors we assume 
plane wave solutions for (j)A — (j^A + 4'^ of the form 

0+ = Ae''^+''+Be-"'+'', 0^ = Ce'"-^ + L'e-*"-^. (28) 

Then Eq. ( [22b| gives (/^ = [-iky±a±]/[e' ~ S\.) 

0± = /±yle'"±^ + /±Se-*"±^, (29) 

with A, B replaced by C, Z?, respectively, if the lower up- 
per — sign is used in and /. Further, Eq. (22a) gives 

= /i±Ae*"±^ + /i±Be-'"±^, (30) 



with/i± = {{e' - if -kl - oij^^\l{t' {e' - Substituting 



Ab' in Eq. (|22d) gives 
9% 



^g^h'^Ae"'^'' +g^h'^Be 
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(31) 



where — [iky±a±]/[e' + 6]; the upper ± signs in / and 
g correspond to the subscripts of a and the lower ones to 
those in front of a. The eigenstates are 
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is a normalization constant, such that each state car- 
ries a unit current, and is given by 



(33) 



2Wa±{t'{e' + (5) + (e' - (5)2 - fcg - a^) ' 

^'^')^ can be 



The solution = (^I^^, ^>b, *b 
rewritten in the matrix form 
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and 
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^ — ta — x 



(36) 



The columns of the matrix product QM are the (unnor- 
malized) eigenstates of our system. 
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